Beyond linear response spectroscopy of ultracold Fermi gases 
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We study RF-spectroscopy of ultracold Fermi gas by going beyond the linear response in the 
field-matter interaction. Higher order perturbation theory allows virtual processes and energy con- 
servation beyond the single particle level. We formulate an effective higher order theory which 
agrees quantitatively with experiments on the pairing gap, and is consistent with the absence of the 
mean-field shift in the spin-flip experiment. 



A single atom interacting with a coherent field displays 
various coherent phenomena, such as Rabi oscillations, 
which can be theoretically described by the exact non- 
perturbative solution of the quantum time evolution of 
the system In many other contexts, however, it is 
sufficient to treat the field-matter interaction within the 
linear response, i.e. to first order in perturbation theory; 
this allows one, for instance, to calculate the response of 
matter in a complicated many-body state 0. In ultra- 
cold Fermi gases of alkali atoms, the interaction between 
atoms and fields may be extremely coherent while, at 
the same time, the atoms can be strongly interacting or 
in a non-trivial many-body state such as superfiuid of 
Cooper pairs. This is likely to lead to new phenomena 
in the field-matter interactions. The first indication of 
this was the absence of mean-field shifts observed (3| in 
the spin-flip experiment using RF fields [4,3]. The ex- 
perimental results were in contradiction with predictions 
given by linear response. However, they were elegantly 
explained 0, ll| as a collective, coherent rotation of the 
Bloch vector within a framework that can be applied to 
Hartree-type mean-field interactions. It has now become 
topical to ask what happens in a similar situation but 
for a more complicated many-body state, such as the 
recently observed Fermi superfiuids p . RF-spectroscouy 
was used for detecting the pairing gap in such systems [3 . 
The results were qualitatively in good agreement |5l^llo| 
with the lowest order perturbation theory. 

In this letter we combine higher order perturbative ap- 
proach of the field-atom interaction with the idea of a col- 
lective coherent rotation of all atoms 0| and apply it to 
the scheme used in the experiment 0- The key point is 
that only total energy conservation is demanded instead 
of energy conservation in each single-atom process. This 
approach allows the transfer of large numbers of atoms, 
as observed in the experiment ^ and is consistent with 
the absence of mean-field shifts in the spin-flip experi- 
ment 1^. It is a step towards high precision description 



of the RF-spectroscopy of superfiuid Fermi gases, and it 
highlights that ultracold Fermi gases may display new 
phenomena in field-matter interactions. 

We study the RF-spectroscopy of a trapped Fermi gas, 
as described in Fig. ^ using perturbative expansion of 
the field-atom coupling. The interaction is described in 
the rotating wave approximation, for a spatially constant 
field with the Rabi frequency 17, by 



Hi = n j dv [e^'^**^(r)«'g(r) + h.c] 



(1) 



where lo is the detuning of the RF-field and the field op- 
erators are expanded in the eigenstates of the harmonic 
trap potential ^'<T(r) = Y.nim.i'niyni^)'^ni7na- In the fol- 
lowing, the term "momentum conservation" should be 
understood, in case of a nonuniform g the conser- 
vation of the trap quantum number. 

The lowest order perturbation theory in Hi using the 
exact trap states gives the transfer probability Pi = 
J dr J dr' P{i,r,r') for atom i = {ni,li,mi). Here 

P(i,r,r')=^ fdn Tdiie'-^*^-^^' Vvl/,,,.(r,r') 
^ -^0 Jo 

X (cig (n ) Cjg 1 ) ) (Cjv , (n )c], ^ 1 ) ) , 

(2) 

where the function ^'y^' (r, r') = tpi{v)\pji{v)tpj{v')ipji{v') 
describes the overlap of atom wavefunctions. The expec- 
tation value for the number of transferred atoms in the 
first-order perturbation theory is obtained by summing 
Eq.(|2l over all atom indices i = {ni,li,mi). In addi- 
tion, the integration over the times ri and ti gives the 
energy conservation condition for a single-atom process. 
We consider the term of the order N in the perturbation 
expansion, (c.f. Fig. |2Jl 
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Figure 1: The RF-field breaks a pair of atoms in internal 
states g' and g by driving a transition to a third state. At 
the single-atom level, the energy cost of creating an excitation 
in the superfluid, 5Ek, has to match the RF-field detuning ui, 
causing a shift of the spectral peak. Higher order perturbation 
theory allows processes where energy is conserved at the many 
particle level, e.g. resonant transfer of the two atoms above 
even when they are individually out of resonance. 
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Figure 2: In the sequential transfer (solid arrows) of A'^ = 4 
atoms, the (time)order in which the atoms are transferred 
can be chosen in A''! different ways. An A^-atom transfer (one 
example shown by dashed arrows) can be done in A^!^ different 
ways. The process is called virtual if an auxiliary vertical line 
(example shown as dotted line) cuts two or more arrows [Tll |. 



We take this to be the expectation value for A^ atoms 
being transferred in a single process. Thus, we neglect 
oscillations of the atoms back and forth corresponding to 
terms of higher order in Hi but with the same final num- 
ber of transferred atoms. In principle, these oscillations 
contribute to the linewidth of the field-atom interaction. 

For weak excitations (linear response), the single- 
particle energy conservation gives the transfer probability 
{5Ek-u))'^+T'^ atom with energy cost SE}^, where lo 

is the field detuning and F is the linewidth. For strong 
single-particle excitations, i.e. Rabi oscillations, this be- 
comes 1 thereby the field-atom coupling Vt 

acts as an effective linewidth. Large effective linewidth 
(F or 17) allows to transfer a large number of atoms non- 
resonantly also in single-particle processes. However, 
the reported linewidth in is small, ~ 100 Hz, which 
could also be deduced from the sharpness of the narrow- 
est features of the spectra, giving a limit to homogenous 
broadening (other features show inhomegenous broaden- 
ing due to pairing and the trap potential). The effective 
linewidth needed for the transfer of the observed num- 
ber of atoms can be estimated, using the calculated dis- 
persion 5Ek, to be ~ 10 kHz. Since this is significantly 
larger than 100 Hz, considerations beyond single-particle 



processes are well motivated. 

Assume an effective pulse length Tp (either the true 
pulse length or the coherence time of the interaction) . If 
the number of transferred atoms A^ is large, the transfer 
processes occur at roughly even intervals. This means 
that each intermediate virtual state lasts for roughly the 
same time tintor — For pulse lengths of = Is, 

and the number of transferred atoms of the order N — 
10^ the corresponding energy uncertainty for any single 
virtual state is /S.E/h > = Kf' Hz (repeated trans- 
fer of smaller A^ within coherence-limited Tp gives the 
same result). Since the typical Fermi energies are of the 
order of 10 kHz, this allows any atom to be transferred 
regardless of its energy cost. However, the total Nth 
order process should conserve energy, meaning that the 
sum of 5Ek — u! oi all transferred atoms should be zero. 

The effective pulse length Tp enters the picture only 
through the uncertainty relation, by introducing an en- 
ergy cutoff 6c for the maximum energy cost of an atom 
transferred through an intermediate (virtual) state. This 
allows to simplify Eq. Q by fixing the cutoff Sc- Thereby 
the actual duration of the virtual state fjnter becomes ir- 
relevant for Winter < J- and we can choose tintcr 0. 
Considering the RF-spectroscopy in the sense of a local 
density approximation, calculating the spectra for each 
position r separately and then integrating over r, we also 
demand that the atoms may not exchange energy in the 
Nth order process over a distance. The probability of 
finding the atoms at the same position is given by the 
overlap of the atom wavefunctions. These approxima- 
tions can be written as 



(Neit)) 



where 
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Hf{t) = C j dr [ne'^'-^l{r,t)-9g{r,t) + h.c.] 



N 



(5) 



and the normalisation constant C = V'^~^. The inte- 
gration over the position r gives the total momentum 
conservation. Having only two time variables in Eq.l^J 
is related to the simultaneous transfer approximation 
(^intcr ~^ 0), and the single position variable in Eq.© 
corresponds to the requirement of the wavefunction over- 
lap. For fermions, the operator Hf as expressed above 
vanishes for A^ > 1 because of the Pauli exclusion princi- 
ple. However, restricting oneself to a subset where each 
one-atom transfer preserves momentum separately, one 
obtains a nonvanishing contribution even for fermions. 
That is, we replace 

^l{r,t)^gir,t) ^ lV'nim(r)|'cV^(i)c„i™g(t). (6) 
nlm 

The process described here is therefore only possible in 
presence of momentum conservation in the field-matter 
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interaction, as is indeed the case in the considered sys- 
tem. 

Our Nth order theory is now simply a linear response 
theory for the defined coherent N-atom coupling operator 
Hf (t). One obtains using the Wick's theorem (when no 
atoms in state e initially) 

Jo Jo 

(7) 

where the preiactor G = y^2jv — i the summation 
is over all A^-atom sets {i} and the probability that this 
set of atoms {«} is transferred is 

^«=E n *.,(r,r')(4(Ti)c,,(ti))(c,e(ri)4(ii)). 

_ (8) 
Here {j} denotes all permutations of the set {i}. Note 

that even though each one-atom process preserves mo- 
mentum, the Nth order process still contains off-diagonal 
Green's functions {clijngCn'i'm'g) ■ 

We approximate the sum over all permutations {j} by 
N independent sums, thus replacing the N\ terms by 
terms. This approximation contains all the correct terms 
but it also introduces additional terms. However, the di- 
agonal Green's functions dominate over the off-diagonal 
ones, and these additional terms give only a small contri- 
bution. In addition, the product over the functions 
gives cancellations of the added terms upon r-integration 
since, for these terms, the atom wavefunctions do not nec- 
essarily appear in the square form l^'ij p. One obtains the 
following form for the transfer probability of set {i} 

i<N j 

(9) 

This probability can be formally written at T = as 

p{^} = n EE^^''''^"''"^^'^'*'^^(*'j'^.r'r'), (10) 

i<N j m 

where AE(i, j,m) is the energy difference of the single- 
atom states and F(i, j, m) is a product of Fermi functions 
and Bogoliubov coefficients giving the occupation and 
transfer probabilities, irrespective of the energy conser- 
vation. The internal variable m, describing the quasipar- 
ticle states, comes from the Bogoliubov-deGennes equa- 
tions. Comparison of Eq.© to Eq.jSJ shows that the 
first-order term contains the same spectral weight fac- 
tors, if we choose j = j' in Eq. ||2J). The sum in Eq. Q 
becomes 

E E e"^^<^''^'"'^^^'~''^F{i,j, m, r, r'). (11) 

j m 

Integration over the times ti and ti gives the single atom 
energy conservation. In (|10|l . the exponential functions 



in the product over i can be combined and upon inte- 
gration over times ti and ti in Eq.Q, yield the total 
energy conservation. This acts as a boundary condition 
for the choice of the transferable atoms {(n^, li^rui)} and 
hence for the actual number of transferred atoms N. The 
expectation number of transferred atoms is now the aver- 
age of all combinations {A} that satisfy the total energy 
conservation. The spectral weights, or the probabilities, 
of different combinations are determined by the prod- 
uct of the terms F(i,j, m, r, r') which also appear in the 
first order spectra in Eq. Hll|) . We approximate the ex- 
pectation value by keeping in {A} only the combination 
that transfers the highest number of atoms. We have 
tested this approximation with 40 atoms and found good 
agreement between the approximative scheme and the 
full combinatorial treatment. In calculating the total en- 
ergy conservation, the energy cutoff Sc for the maximum 
energy cost of an transferrable atom acts as an additional 
constraint for the possible combinations. 

We now use the first order spectrum to calculate the 
number of transferred atoms by the following algorithm. 

• Consider a fixed field detuning uj. Set 6min — —oo, 
^max = oo. Let Ni{d) be the expectation number 
of atoms that one can transfer using the first-order 
perturbation theory and by giving energy 6. 

• Integrate dS Ni{S){S — lu) =: 
AE{drmn,6max,i^)- This givcs the total en- 
ergy change if all the atoms with energy change of 
S e [Srmn,Srnax] are transferred. 

• If AE{Srmn,Smax,i^) IS positive, decrcase Smax or 
if it is negative, increase 6min and go back to the 
previous step. Otherwise, proceed to the next step. 

• The maximum number of atoms that one can trans- 
fer is NNiuj) J^2T d^Nii^)- 

One can map the whole spectrum by using the algorithm 
above for several detunings uj. For a nonuniform gas, the 
process should be repeated for all positions r and r' and 
integrated to give the total number of transferred atoms. 

For a noninteracting gas, the first-order perturbation 
theory and the current higher order theory give identical 
results since the conservation of the total energy equals 
the conservation of energy of a single one-atom process. 
For atoms interacting only via the Hartree field, the two 
theories differ by the position of the spectral peak. As an 
example we consider the spin-flip experiment of Ref. j3| , 
with Ng {Ng' ) atoms in the g-state (g'-state) and mutual 
interaction energy Vgg>. In the first-order perturbation 
theory, the RF-field needs to give the energy correspond- 
ing to a single atom excitation, i.e. AE = Vgg' {Ng — Ng'). 
The higher order theory requires only the energy cor- 
responding to the complete iV-atom process. Since the 
mutual interaction energy of the two fermion species g 



4 



and g' is unchanged if the spins of all atoms are flipped, 
the total energy change of the whole process is zero [l^ ■ 
This is consistent with the absence of the mean-field shift 
observed in the spin-flip experiment in Ref. |^. 

To describe the superfluid Fermi gas in a harmonic 
trap, we apply the mean-field resonance superfluidity 
approach used in Ref. Q but restricting oneself to the 
one-channel model in order to describe broad Feshbach- 
resonances 0, We calculate self-consistently the 

position dependent excitation gap, density distribution 
and the chemical potential for 2670 atoms using back- 
ground interaction strength U = 2.itkvQr^^^ and cutoff 
Ec = 161.5 fiwo, where fiwo, Tosc are the oscillator en- 
ergy and the length. The position dependent excitation 
gap is 0.83 Ey at the center of the trap corresponding 
roughly to Li^ atoms with peak density of 10^'^1/cm^ 
at the magnetic field oi B k, 834 G, as obtained from 
the renormalised theory for a uniform gas Q . The small 
number of atoms used above gives too narrow density and 
gap profiles, and one expects a larger excitation gap close 
to the edges of the trap for numbers of atoms used in the 
experiments. On the other hand, for pulses longer than 
the coherence time of the atoms, the system can relax 
and lower the excitation gap during the pulse (not con- 
sidered here). Thus, one expects the two approximations 
to partly cancel each other. 

Fig. O shows the RF-spectra at two temperatures T = 
0.0 T-p and 0.08 Tp using a 7r/2-pulse. The peak is shifted 
due to the pairing gap and, most notably, another peak 
at the zero detuning exists at a finite temperature, origi- 
nating from contributions from the edges of the trap [Sj . 
This behaviour and the locations of the two peaks at 
finite temperatures are in good quantitative agreement 
with the experiment Note that, in principle, first- 
order perturbation theory gives a maximum transfer ra- 
tio of only ^ 5% of the atoms, because the single particle 
excitation energies 5Ek have a strong momentum (fc) de- 
pendence and therefore only a small fraction of the atoms 
can match the energy of a narrow lincwidth field. The 
higher order theory is not limited in this way and al- 
lows to describe the experiments with large numbers of 
transferred atoms, such as 50% in 0; c.f. Fig. |21 

In conclusion, the connection between the first-order 
and higher order terms in perturbation theory was used 
to define an effective higher order theory for the RF- 
spectroscopy of Fermi gases. It was shown to be consis- 
tent with the temperature dependence of the pairing gap 
spectra and the coherence phenomenon observed in 
Ref. [ij. This approach could be applied also to other 
types of field-matter interactions and spectroscopies in 
ultracold gases. It should be useful whenever the single 
particle excitation spectrum is non-trivial, e.g. BCS-type 
states or exotic states in optical lattices, and the process 
is coherent. Due to the special nature of field-matter in- 
teractions discussed here, the understanding of phenom- 
ena such as Josephson oscillations in Fermi gases — when 
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Figure 3: RF-spectra in the higher order theory using param- 
eters given in the text. At finite temperatures, an additional 
peak appears at the zero detuning. For comparison, the spec- 
trum from the linear response theory is shown in solid line. 



the effective tunneling is realized by fields — will be 
non-trivial and may show phenomena not existent in the 
corresponding solid state system. 
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